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对 称 和 矩阵 反问 题 的 总 体 最 小 二 乘 解 ” 
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摘 要 : 最 小 二 乘法 是 近年 来 求解 对 称 和 矩阵 反问 题 的 一 种 常用 方法 ， 但 因 系数 矩阵 常常 存在 误差 ， 方 法 本 
身 具 有 很 大 的 局 限 性 。 鉴 于 此 ， 本 文 提出 并 讨论 了 对 称 矩 阵 反 问题 的 总 体 最 小 二 乘 解 ， 给 出 了 解 
的 一 般 表 达 式 ;证 明了 最 佳 逼 近 问 题解 的 存在 唯一 性 ， 给 出 了 其 具体 表达 式 及 数值 算法 ， 并 将 数 
值 结 果 应 用 于 求解 对 称 和 矩阵 反问 题 。 
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用 Rm 表示 n xm 实 和 矩阵 的 全 体 ; R” = R”; SR” 表示 n 阶 实 对 称 矩 阵 的 全 体 ;O Rx” 
表示 n 阶 正 交 和 矩阵 的 全 体 ， 了 表示 上 阶 单位 阵 ，A+ 表示 和 矩阵 4 的 Moore-Penrose 广 义 逆 ，ImA 
表示 矩阵 4 的 值 域 ， 即 In4 = {hz |x € R"}， 其 中 4A e R"*"; ac(4) 表 示 和 矩阵 4 的 特征 值 的 
全 体 ; tr(4) RIRE ABI; RA 表示 由 4 的 列 向 量 组 成 的 线性 组 合 的 全 体 ， 即 矩阵 4 的 
值 域 。diag( 和 1, 和 2, …，》n) 表 示 对 角 元 素 为 Ai, Xa, … ,和 的 n 阶 对 角 和 矩阵 ; E Rm EE 
义 矩 阵 4 与 妃 的 内 积 为 (4,B) = tr(87A)， 则 由 此 内 积 导 出 的 范 数 上 |4l| = VA, A) 是 矩阵 
的 Frobenius 范 数 ， 并 且 R^*"^ 构成 一 个 完备 的 内 积 空间 。 

文 [1] 讨 论 了 对 称 和 矩阵 反问 题 的 最 小 二 乘 解 。 文 [2] 从 数值 分 析 角 度 研 究 了 总 体 最 小 二 乘 问 
题 。 文 [3,4 和 通过 分 析 比 较 最 小 二 乘 与 总 体 最 小 二 乘 问题 ， 理 论 上 说 明了 总 体 最 小 二 乘 是 比 最 小 
二 乘 更 合适 的 拟 合 方法 。 

本 文 研究 如 下 问题 : 

问题 I AEX, Be m", 求 4e SR"x"， 使 得 AX = B, Kp X, BWE 


min||[X — X, B — B]|, s.t. rank(X)= rank ( ) 
B 
HSR"*"—(A|Ae R?*", AT = A). 
BRI 给 定 4e Rnxn， 求 4e S4， 使 得 
JĀ- Al = inf LÀ - Al, a) 
其 中 54 是 问题 1 的 解 集合 。 
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AB e HT BBLUG EMITSA AR ME, EAIB T MEI — RROA SS HERT AARE 
逼近 解 的 存在 唯一 性 ， 给 出 了 其 表达 式 及 求解 算法 ， 最 后 将 数值 结果 应 用 于 求解 对 称 矩 阵 反问 
题 。 

2 问题 工 的 解 


首先 引入 几 个 概念 。 
定义 2.16 设 M 是 R"? 的 一 个 子 空间 ，M C R"*, 称 M 是 也 - 中 性 的 ， 如 果 


[ry] = (Hz,y) =y Hz =0, Vr,yeM. 


5EX,.2.20501. xpEmxnsSBEEX. nre. id 
In 
G(X) 2 Im e rio 
X 


我 们 把 n 维 子 空 间 G(X) 称 作 关 的 图 。 进 一 步 ， 对 于 Rm"t* 中 的 n 维 子 空间 S$S， 若 S 是 某 
Am x n ERRER, MAS ARFER. 

下 面 考 虑 矩阵 方程 CY zc D(C, D € R"*"), © Hir A iA R B A aD R n] BL 2 
RY e SR”, 使 GY = 万， 其 中 心 ,万 满足 


min||C ~ Ĉĉ, D ~ Dj|, st. R(D)c R(Ó) 


AH Lagrange RTIA, X [5] EB] TEEDE CY 万 带 对 称 约 束 的 总 体 最 小 二 乘 问题 可 
归结 为 求 Riccati 和 矩阵 方程 XSX -- RX - XR—- S — 0 的 解 ， 其 中 R= CTC- DTD, S = 
CT D — DTC。 该 方程 的 解 与 矩阵 


的 不 变 子 空间 密切 相关 。 记 


TÆRNE, HKT = -TK， 则 存在 非 奇异 实 和 矩阵 W 使 得 
W-'!TW = diag[Mx, 一 Xe ，》 —A1, 0, +++, 0], 


其 中 X > M4_1 >.…> M > OË T HERES. 
引 理 2.1 回 设 C DERmxn， 令 


R S 
gm ; 
HKiRR-CTC-DTD, S=CTD+DTC， 且 非 奇 异 实 矩阵 信使 得 


W-'TW 一 diag[àx, —Àk, "T y Ài, 一 入 1， 0, STAR 0J, 
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EP Ak > Xx-_1 > … > 和 1 > 0 是 TT 的 正 特征 值 ， 则 存在 nw 维 的 、 工 不 变 的 、K- 中 性 的 子 空 
间 ,So， 使 得 o(T | So) = {Xk,… ,入 1,0}; 或 如 果 k =n, Hol | So) = Qu A} 

3|z82.20 车 引 理 2.1 中 So 是 图 子 空 间 ， 则 由 它 可 确定 关于 [C, D] 的 带 对 称 约束 的 总 体 最 
小 二 乘 解 ， 特 别 地 ， 当 =n 时 ， 带 对 称 约束 的 总 体 最 小 二 乘 解 是 唯一 的 。 

5178 2.35] 带 对 称 约束 的 总 体 最 小 二 乘 问题 中 修正 矩阵 C 满 足 C = (C+DY)(I+Y?)-!， 
其 中 Y 是 对 称 总 体 最 小 二 乘 解 。 

利用 带 对 称 约束 总 体 最 小 二 乘 问题 的 结果 ， 可 得 对 称 矩 阵 反问 题 的 总 体 最 小 二 乘 解 。 

定理 2.1 设 X, BeEeRnxm， 令 


R S 
r={ ) R=XTX - BTB, S- XT B4 BT x, (2) 
S -R 
HAgdedE S 5e KEE W 使 得 


W-'TW = diagD, —A, ++, Ài, 735 0, 77: , 0), 


HF Ak > àp- > o > A > 0 是 荆 的 正 特征 值 ， 由 W 的 奇数 列 组 成 的 子 阵 ( ) 车 Y € 
Re”? 非 奇 异 ， 则 4 = ZY-! 为 问题 1 的 一 个 解 ， 且 问题 1 解 的 一 般 表 达 式 为 
A= Ao+(IT— XR+)M(I— XX*) vVMeSR"", (3) 


其 中 ho = ZY-!, X = (I+ A0?)-1(X + AoB). 
证 明 设 问 题 I 的 任 一 解 4e SR2xn， 则 有 74 = BT7， 其 中 闷 , 户 满足 


min [X7 一 XT BT- TII, st. R(BT)c R(XT), 


BI 4 可 看 作 关于 [X7, BT] 的 对 称 总 体 最 小 二 乘 解 。 
在 引 理 2.1 中 令 C = XT, D— BT, 得 


R= XTX -BB, S$S= XiB+BTX. 


由 引 理 2.1 知 ,从 阵 ( ) 的 列 张 成 的 子 空间 即 为 满足 引 理 2.1 条 件 的 子 空间 So 


XA So 是 图 子 空间 的 充 要 条 件 是 Y-1 存在 ， 所 以 由 引 理 2.2 知 4 = ZY-1 为 问题 1 的 解 。 

5432823, PF HmUEBA-(I-XX*)Mu-XX*)RJIjg AX =0(A € SR"x") 的 
通 解 。 

事实 上 ， 对 任意 的 M e SR”, A = (I- XX*)M(I - X*)&ZjfRAX = 0(4 € 
SR”) 的 解 ， 而 对 AX = 0 的 任意 一 个 解 4e SR, 令 M = A， 有 

üu-XX*)àu-XX*)-ü-XX*)A- A- XX* AS A- (X*) (S) A- A, 

所 以 A= (— XX*)M(I - R+) 是 方程 AX =0(A € SR"x") 的 通 解 。 

于 是 ， 问 题 I 解 的 一 般 表 达 式 为 

A-Ao-(I- XX*)Mü -XX*) YM ESRnxn， 


其 中 Ao = ZY-1, X = (I + Ao?) (X + AB) 
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3 问题 开 的 解 
313.16] Ec R”, WpEdn ERRER H E 
E+ET 
|e- —7—] s e - en. 


定理 3.1 对 给 定 的 4 cR, ULLA RE GESE2.14RTEI X, B € R^*", HEIFER 
解 Ae Sa bn X neg 5e (B 2 CON 


其 中 

U-(U, Us) Ee OR"™Y", V=(W V) eog*n, 

Ue R, VW eR"*", r—rank(X), 

: A A 
UT(Ã- AU =| "19. ane. 
421 A22 
DJE TI RO f A 可 表示 为 à 
A= Ao + 5Ux (Azo + ÃZ UT. (4) 
证 明 显然 


Sa = {4o + (I —- KK*)Mü - XX*)| Ao = ZY , VM e SR^*"Y 


是 完备 内 积 空间 Rnxn pa AME, HUEHSGEUO GE RRUI 知 存在 唯一 的 À c S4， 使 式 (1) 成 
AE 


下 面 证 明 问 题 开 的 解 4 可 由 (4) 式 表示 。 对 4e S4, id 


Mi Miz 


UT MU = 
M21 M22 


| , Mii € R'™", 


则 


14-Al? = |[A — Ao — (U2U27) M (UaU27)|" 


-Jaw (7 Ai ee ( E )r 


= [v7 (à - Av - m )r 
0 M 


= |Aailf? + HA2ll? + lall? + l| Ã22 — Mozl?. 
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由 引 理 3.1 可 知 ， 当 i 
A22 + Á 

Mam 2z (5) 


tf, JĀ- 4 达到 最 小 值 ， 由 (3) 和 (5) 式 ， 可 得 
4=4o 十 3Ua( Az + 422)D27. 


算法 3.1 求 问 题 开 的 解 的 步 又 如 下 : 
步 1: 按 (2) 式 计算 T，; 


2: 由 定理 2.1 计算 问题 I 的 一 个 解 = ZY 71, BSEC ( 7), JERRY C 是 否 存在， 
步 3， 求 总 ， 并 对 其 进行 奇异 值 分 解 ; 
54s GRUB (4) Sk À. 

4 ”数值 算 例 


下 面 通 过 例子 说 明 前 述 理 论 和 算法 的 应 用 。 
Bj 设 n = 4, m = 3， 给 定 如 下 的 X, Be R**m，A e Rnxn， 求 问题 开 的 解 。 


2.5 3 0 3 1.5 0 3 —2 0 2 
0 1 1 0 0 L5 - 2 1 0 -1 
X= , B= , A4= 
-] 0 05 -1 0 0 1 2 1 -3 
15 2 0 3.5 1 0 -1 1 2 2 
计算 
R S 
T- 
S -R 
其 中 
4 3 0.5 -235 24 15 —5.5 19.25 
3 一 0.25 0.5 2 1.5 3 0.75 1 
R= ; S= 
0.5 0.5 0.25 2 —5.5 0.75 2 一 5 
一 2.25 2 2 一 了 19.25 1 一 5 14.5 
] Y 
计算 ( ) xm 
Z 


0.1038 — 0.1724 0.1688 0.5806 

0.1205 — —0.0552 0.5822 0.0812 

0.1839 . 0.3568 0.4284  —0.1108 
—0.2323 0.0087 0.0104 0.3340 
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0.3573 
一 0.1417 

0.6450 

0.5685 

易 知 Y-! 存在， 于 是 计算 

0.0033 
一 0.8021 

1.1531 

2.0362 


计算 广 ， 并 对 其 进行 奇异 值 分 解 
3.0053 — 2.4692 
: 0.1400 — 0.8430 

X 


—1.0436 0.2744 
1.8263 1.2674 


—0.8510 . 0.0347 
—0.1496 0.7466 


li 


0.1365 . 0.6557 


0.4740 
—0.4588 
0.4457 
—0.4539 


—0.8021 
1.7844 
—0.6386 
0.6716 


0.1909 
1.0500 
0.5502 
0.0091 


—0.2555 
0.6476 
—0.7163 


—0.4846  —0.1066 0.0470 


4.5746 0 0 
0 1.5330 0 
0 0 0.3607 
0 0 0 
利用 (4) 式 计算 4 得 
1.3729  —0.7218 
Á- —0.7218 1.7892 
0.5671 — —0.6731 
—0.5585 0.5191 


0.0486 

0.6369 

0.0772 
—0.1870 


1.1531 
—0.6386 
0.6442 
—2.0827 


—0.4576 
—0.0269 
0.1958 
0.8669 


—0.7882 
—0.6130 
—0.0544 


0.5671 
—0.6731 
0.8950 
—0.9725 


0.4891 
—0.0257 
—0.1496 

0.5198 


2.0362 
0.6716 
—2.0827 
—2.8378 


—0.4372 
0.4957 
0.7504 


—0.5585 
0.5191 
—0.9725 
2.0777 


0.4330 
—0.6153 
0.6587 
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通过 求解 得 : 上 ArLs 一 Á|| = 6.5317. ||Àrs — Á|| = 10.8546, X&'P Áris. Ars 为 分 别 运用 
AS PE R^] REMER RERET. 
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The Total Least Squares Solution for the Inverse Problems of 
Symmetric Matrices 
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Abstract: Least squares have been widely used in the inverse problems of symmetric matrices. How- 
ever, errors always occur in the relevant coefficient matrix. This makes the approach limited. In order 
to overcome this shortcoming, the total least squares solution of inverse problems of symmetric matrices 
are proposed. The general form of the solution is given. The existence and expression of the optimal 
approximation solution are presented, and a numerical algorithm is derived. These results are finally 
applied to solve the inverse problem of symmetric matrices. 

Keywords: symmetric matrix; inverse problem; total least squares solution; singular value decompo- 
sition; Riccati equation 


